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ABSTRACT 


This  is  a final  report  of  SOFAR  (Sound  Fixing  And  Ranging)  float 
trajectory  data  from  MODE  (Mid-Ocean  Dynamics  Experiment)  for  the 
period  September  1972  through  June  1976.  It  begins  with  a brief 
historical  review  and  includes  a summary  of  the  steps  involved  in 
processing  the  data  and  computing  the  trajectories.  The  data  is 
presented  in  three  ways:  1.  Time  annotated  trajectories  on  a 
Mercator  projection;  2.  Stick  plots  of  float  velocity  as  a 
function  of  time;  and  3.  A sequence  of  charts,  each  spanning  ten 
days,  showing  the  spatial  relationship  and  relative  motion  of  all 
floats.  There  is  also  a table  of  launch,  recovery,  depth,  number 
of  days  tracked,  average  velocity,  speed,  and  kinetic  energy  for 
each  trajectory. 


INTRODUCTION 


The  concept  of  using  the  SOFAR  channel  to  acoustically  locate 
neutrally  buoyant  floats  at  great  distances  was  first  proposed  by 
Stommel  (1949)  and  later  advocated  by  him  in  a letter  to  Deep-Sea 
Research  (1955).  His  original  idea  was  based  on  the  use  of  SOFAR 
DEVICES  (a  euphemism  for  explosive)  that  would  be  released  from  a 
carrier  float.  For  understandable  reasons  this  technology  was 
never  developed. 

In  1965,  a program  of  development  was  initiated  by  Stommel,  Webb, 
Tucker,  and  Volkmann  to  evaluate  the  feasibility  of  long-range, 
very  low  frequency  pulsed  CW  transmission  in  the  ocean  (Tucker  and 
Webb,  1970).  Encouraging  results  stimulated  the  development  of  a 
neutrally  buoyant  SOFAR  float  (conceptually  identical  to  the 
Swallow  float),  which  was  first  used  successfully  in  1969  in  the 
Bermuda  Triangle  (Rossby  and  Webb,  1969;  1971). 

Not  long  after  this,  Stommel  announced  plans  to  organize  a synoptic 
study  of  ocean  currents  and  their  variability  in  the  western  North 
Atlantic.  This  program  later  came  to  be  known  as  MODE.  A proposal 
was  submitted  to  the  International  Decade  of  Ocean  Exploration 
Office  at  the  National  Science  Foundation  by  Webb  and  Rossby  to 
develop  a new  and  less  expensive  type  of  SOFAR  float.  From  a 
technological  point  of  view  this  was  a high  risk  program,  since 
the  new  float  design  had  not  been  tested  and  the  timetable  had 
already  been  prescribed  by  the  field  program  as  a whole.  Never- 
theless, the  proposal  was  funded  and  the  float  design  proved  to 
be  successful.  Twenty  floats  were  built  for  MODE.  The  details 
of  the  field  program,  its  operation,  and  some  of  the  data  that 
were  collected  have  been  reported  in  the  literature  (Rossby,  Voorhis 
and  Webb,  1975;  Freeland,  Rhines  and  Rossby,  1975)  and  will  not  be 
repeated  here. 
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This  document  contains  plots  of  all  float  trajectory  data  that  have 
been  collected  between  September  1972  and  June  1976.  In  addition 
to  the  trajectory  plots,  the  velocity  of  each  float  along  its  path 
has  been  computed  and  plotted  as  a time  series  stick  plot.  A series 
of  ten  day  charts  for  the  period  April  1973  through  January  1974, 
showing  the  relative  motion  of  all  floats  as  a function  of  time,  is 
also  included.  For  reference,  elementary  statistics  (i.e.,  number 
of  days  tracked,  average  speed,  velocity,  and  kinetic  energy)  are 
shown  in  Table  1,  but  the  reader  is  cautioned  that  the  averages  may 
carry  little  significance  if  the  trajectory  spans  different  geo- 
graphical regions  which  have  different  energy  levels  and  dynamics. 

In  addition,  the  depth  of  each  float  after  launch  is  given.  Due 
to  compressional  creep  of  the  main  cylinder,  floats  sink  at  a rate 
of  0.85  meters/day  (Rossby,  Voorhis  and  Webb,  1975). 


SIGNAL  TRANSMISSION  AND  RECEPTION 


The  MODE  experiment  was  centered  at  28°N,  69°40‘W.  All  floats  were 
launched  in  this  area  and  ballasted  for  1500  meters  depth  (slightly 
below  the  sound  channel  axis).  (Figure  1 shows  a float  being 
launched  from  the  R/V  RESEARCHER).  The  acoustical  signals  trans- 
mitted by  each  float  radiate  horizontally  through  the  ocean  at  a 
speed  of  about  1492  meters/second.  Their  time  of  arrival  (TOA)  at 
each  of  several  stations  yields  distances  from  the  float,  provided 
that  the  transmission  time  is  known.  The  primary  requirement  for 
locating  a float  is  to  determine  the  TOA  correctly.  This  is  done 
by  using  receivers  with  a very  sharp  bandpass  filter  carefully 
tuned  to  the  acoustic  carrier  frequency  and  by  limiting  the  search 
for  the  signal  to  a narrow  window  in  time  centered  near  the  expected 
TOA. 

Four  stations  have  been  used  in  MODE  to  track  the  floats.  They  are 
located  at  Bermuda,  Eleuthera  (Bahamas),  Puerto  Rico,  and  Grand 
Turk.  Except  for  the  Bermuda  station,  all  are  equipped  with  digital 
systems  that  continuously  store  the  data  on  computer  tape.  All 
stations  are  equipped  with  backup  systems  which  record  the  signals 
graphically  on  a drum  recorder.  The  receiver  at  Grand  Turk  is  shown 
in  Figure  2. 

Three  acoustic  carrier  frequencies  are  employed  with  up  to  seven 
floats  operating  on  each  frequency.  These  floats  are  identified 
by  their  pulse  repetition  rates  (all  of  which  are  close  to  one/ 
minute).  In  order  to  make  the  transmission  schedule  repeat  each 
day,  they  are  arranged  to  transmit  exactly  1437,  1438,  1439,  1440 
(one/minute),  1441,  1442,  or  1443  pulses/day.  Thus,  on  a drum 
revolving  one  turn/minute,  each  float  has  a different  slope  making 
it  easy  to  recognize.  A sample  drum  recording  with  each  float  TOA 
labeled  is  shown  in  Figure  3.  The  float,  channel,  and  period 
information  is  summarized  in  Table  2. 
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Figure  1 A SOFAR  float  is  being  launched  from  the  R/V 
floats  are  stored  on  deck  for  future  launches. 


Three  other 


Figure  2. 

Above: 

Listening  station  at  Grand  Turk. 
Electronic  equipment  is  housed  in 
the  small  center  shack. 


Left: 

Grand  Turk  float  signal  receiver 
and  recording  equipment. 


Photographs  by  David  J.  Dorson 


Channel-Period  Table 


tdcn  receiver  consists  ot  a permanently  mstaiieci  riyurupnunc,  at 
the  depth  of  the  SOFAR  channel,  connected  by  a shore  cable  to  the 
amplifiers,  filters,  and  recording  equipment.  The  preamplifiers 
are  very  high  gain  linear  systems  with  high-pass  filters  near  the 
input  to  suppress  any  harmonics  (due  to  non-linearities)  that 
might  arise  from  the  energetic  low-frequency  shipping  noises.  The 
output  is  passed  to  high-Q  filters,  which  bandpass  each  of  the 
three  carrier  frequencies.  The  3 dB  and  equivalent  noise  band- 
widths  are  0.7  and  1.2  Hz  respectively.  After  fullwave  rectifi- 
cation, the  outputs  are  multiplexed  into  an  analogue  to  digital 
converter,  which  samples  and  ■'•ecords  each  channel  five  times  per 
second  as  an  11-bit  word  on  a computer  compatible  seven-track, 

556  BPI  tape. 


DATA  PROCESSING  PROCEDURES  AND  PROGRAMS 


The  data  processing  was  done  at  the  Yale  Computer  Center  on  the 
IBM  1401  and  7040-7094  DCS  system  and,  later,  on  the  IBM  370-158. 
The  float  tracks  were  plotted  at  Yale  and  the  stick  plots  at  the 
University  of  Rhode  Island  on  an  HP2000  multiprogramming  system 
using  a Versatec  plotter.  The  three  major  computer  programs  used 
are:  GRIND-tape  preprocessor,  CRUNCH-T0A  finder,  and  TRACK-float 
position  calculator.  The  GRIND  program  was  written  in  ASSEMBLER, 
and  the  others  were  written  in  FORTRAN  IV. 

The  GRIND  program  reads  the  seven-track,  556  BPI  tapes  mailed  from 
the  stations,  reformats  the  data,  and  writes  a nine-track,  1600  BPI 
tape.  Data  is  checked  for  correct  number  of  channels,  correct 
order  of  channels,  missing  data,  data  off  scale,  and  irregular 
station  clock  increments.  Any  records  having  errors  are  flagged 
on  the  output  tape  and  printout. 

CRUNCH  reads  the  tape  output  of  GRIND,  finds  the  TOAs  for  each 
float,  and  prints  and  punches  the  TOAs.  The  tape  is  read,  and 
any  data  occurring  during  a float  window  is  added  to  the  summed 
window  for  that  float  (see  Figure  4).  At  the  end  of  a four-hour 
period  of  data,  the  TOAs  are  found  by  correlating  each  summed 
window  with  the  expected  signal  shape. 

After  several  weeks  of  TOAs  have  accumulated,  data  is  checked  by 
hand.  Any  datum  that  is  missing,  or  obviously  bad,  is  replaced  by 
taking  data  from  the  backup  drum  recordings  or,  when  possible,  by 
interpolating  data  from  surrounding  good  points. 

The  basic  element  in  the  TRACK  program  is  to  compute  the  great 
circle  distance  between  two  points  on  the  Clark  spheroid  using 
the  Andoyer-Lambert  correction  (Thomas,  1965).  The  procedure  is 
iterative.  For  circular  navigation  (also  known  as  range-range), 
a position  for  the  float  is  assumed  and  the  distances  to  the 
receivers  computed  and  then  compared  to  the  observed  values 


Computer  tape 
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(Distance  = Travel  time  x Speed  of  sound).  The  errors  are  used 
to  estimate  a new  position,  and  the  process  is  repeated  until  the 
sum  of  errors  is  less  than  100  meters.  The  process  of  computing 
a new  position  is  accomplished  by  comparing  the  travel  time 
errors  to  those  created  by  mapping  a perturbation  vector  in  geo- 
graphical coordinates  into  travel  time  perturbations. 

The  same  iterative  procedure  is  also  used  for  hyperbolic  tracking 
except  that  the  quantities  now  being  compared  are  the  arrival  time 
differences  at  pairs  of  stations.  In  summary,  the  advantage  of 
the  iterative  scheme  is  its  simplicity  and  the  ease  with  which  it 
could  be  programmed. 


TRACKING  ACCURACY 


Throughout  MODE  there  were  nearly  50  launches  and  35  recoveries  of 
SOFAR  floats.  The  position  of  each  of  these  operations  is  known 
from  the  shipboard  satellite  navigation  systems.  The  correspond- 
ing position  of  each  float  can  be  determined  using  the  previously 
described  triangulation  methods.  Assuming  the  satellite  fix  is 
without  error,  it  is  interesting  to  see  how  they  compare.  Before 
doing  this,  however,  we  shall  summarize  briefly  several  sources 
of  error  that  may  affect  the  tracking  accuracy.  These  are: 
speed  of  sound  variations,  oscillator  offset  and  drift,  baseline 
magnification,  uncertainty  in  TOA  calculations,  and  receiver  clock 
errors. 

The  horizontal  variations  of  the  sound  velocity  minimum  are  quite 
small  (<  ± .2%),  except  close  tn  the  Caribbean  islands  where  the 
shoaling  of  the  thermocline  and  the  Antarctic  Intermediate  Water 
give  rise  to  a shallower  and  lower  sound  velocity  minimum  (Rossby 
and  Webb,  1969).  The  temporal  variations  are  probably  somewhat 
smaller,  because  the  speed  of  sound  is  a quantity  which  has  been 
averaged  over  several  hundred  kilometers,  spanning  several  extrema 
of  the  spatial  variability.  We  believe  the  temporal  errors  do  not 
exceed  ± .1%  and  probably  occur  on  "eddy",  or  weekly,  time  scales. 
Thus,  throughout  MODE  we  have  assumed  a value  of  1492  meters/ 
second  everywhere.  We  have  no  way  of  independently  monitoring 
the  sound  velocity  variations  from  the  floats  to  each  station. 

Oscillator  offset  and  drift  are  a major  source  of  error,  but  these 
can  be  kept  under  control.  By  hyperbolic  tracking,  the  clock 
drift  can  be  monitored  and  later  corrected  in  the  tracking  program. 
After  offset  and  drift  correction,  the  remaining  clock  errors  are 
usually  < ± 1 second,  and  typically  ± 0.5  seconds,  which  probably 
reflect  speed  of  sound  uncertainties  since  crystal  oscillators 
tend  to  drift  monotonically.  In  fact,  this  error  could  be  re- 
duced to  a negligible  amount  if  the  speed  of  sound  were  indepen- 
dently known  to  sufficient  accuracy. 
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Baseline  magnification  errors  occur  if  a float  is  much  farther 
away  from  the  receivers  than  the  distance  separating  the  receivers, 
or  if  the  float  is  close  to  the  baseline  between  the  receivers  or 
an  extension  of  the  baseline.  Errors  magnify  as  the  angle  of  inter- 
section between  the  equidistant  circles  (or  correspondingly  the 
hyperbolae)  decreases.  This  is  no  problem  for  floats  in  the  MODE 
area,  but  is  troublesome  for  floats  to  the  north  or  west  or  very 
close  to  the  Caribbean  islands.  In  addition,  several  hydrophones 
can  be  shadowed  from  the  float  by  surrounding  islands,  in  which 
case,  there  is  no  way  of  independently  checking  the  calculations. 

There  is  also  uncertainty  in  calculating  the  TOAs.  One  way  to 
obtain  an  estimate  of  this  uncertainty  (or  at  least  a lower  bound) 
is  to  look  at  the  scatter  between  two  signals  which  are  transmitted 
with  a fixed  delay  six  minutes  apart.  For  strong  signals,  the 
standard  deviation  for  this  scatter  is  about  0.25  seconds.  If 
this  variance  is  attributed  equally  to  the  main  and  delay  pulse, 
the  standard  deviation  for  a single  strong  TOA  is  about  0.2  seconds, 
which  is  equivalent  to  the  sampling  rate  of  5 times/second.  As 
the  signal  to  noise  ratio  decreases,  this  jitter  increases  to 
0.6  seconds  and  more.  In  this  case,  it  becomes  necessary  to  edit 
the  TOA  data  sets  before  they  can  be  used  for  tracking.  One 
should  remember,  however,  that  the  above  discussion  does  not 
include  errors  due  to  either  speed  of  sound  variations  caused  by 
internal  waves  or  changes  in  the  signal  shape  on  time  scales 
longer  than  six  minutes.  It  is  our  impression  that  strong  signals 
(no  noise  contamination)  are  relatively  unaffected  by  the  complex 
multipath  arrival  structure  which  is  so  well  documented  for  the 
classical  SOFAR  shot.  This  simply  means  that  most  of  the  "shot" 
energy  is  concentrated  at  the  end  of  the  received  signal.  Thus, 
from  our  point  of  view,  the  impulse  response  of  the  sound  channel 
is  nearly  an  impulse. 

Clock  errors  at  the  receivers  are  negligible,  since  they  were 
frequently  checked  against  WWV. 

In  summary,  for  floats  in  the  MODE  area,  the  uncertainties  are: 

< ± .1%  (^  ± 0.5  sec.)  for  the  speed  of  sound,  ± 0.5  sec.  for 
the  oscillator  offset  and  drift,  a small  multiplicative  factor 
in  the  MODE  area  for  baseline  magnification,  ± 0.2  - 0.6  sec. 
for  TOA  calculations,  and  negligible  for  the  receiver  clock 
errors.  The  cumulative  effect  of  these  is  < ± 1 second. 

Returning  to  the  promised  comparisons,  we  show  in  Figure  5 the 
circular  float  fixes  relative  to  the  ship  fix  for  a number  of 
launches.  These  measurements  are  the  difference  between  the 
float's  position  alongside  the  ship  on  the  surface,  and  the  first 
few  fixes  via  underwater  sound.  These  are  given  for  the  hydro- 
phone pairs  Eleuthera/Grand  Turk  and  Eleuthera/Puerto  Rico  for 
clocks  that  are  adjusted  hyperbol ical ly  immediately  after  launch, 
and  clocks  where  no  adjustment  is  made.  There  is  less  scatter 
to  the  float  fixes  if  the  clock  is  corrected  hyperbol ical ly. 

The  overall  scatter  of  the  float  fixes  is  about  ± 2 km.  Bearing 
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Accuracy  of  Float  Positions  (Km) 


in  mind  that  each  dot  is  a separate  launch  at  a different  time 
and  place,  we  feel  the  result  is  rather  good  and  not  inconsistent 
with  the  above  error  discussion.  The  offset  of  the  float  fixes 
to  the  northeast  relative  to  the  ship  is  due  to  the  assumption 
that  the  speed  of  sound  is  1492  meters/second  everywhere.  Of 
course,  significant  deep  motion  or  shear  across  the  thermocline 
can  introduce  additional  errors. 


CALCULATION  OF  FLOAT  VELOCITIES 


The  velocity  of  a float  along  its  track  is  defined  as  its  rate  of 
change  of  position  as  a function  of  time.  Since  the  individual 
positions  have  uncertainties  on  the  order  of  1 km,  it  is  not 
meaningful  to  compute  velocity  estimates  every  few  hours.  Instead, 
the  time  series  of  position  is  low  pass  filtered  and  subsampled  to 
give  velocities  once  a day. 

The  procedure  is  straightforward:  given  an  array  of  2N+1  positions, 
xi , as  a function  of  time,  ti , representing  either  the  N/S  or  E/W 
coordinate,  we  calculate  the  low  pass  filtered  velocity  by  fitting 
in  a least  square  sense  a parabola,  xi , to  the  array.  That  is,  we 
define  the  error  function 


i=+N 


,2  _ 


= E (Xi-Xi)2  = 


E (xi -a-bt-j -ct -j  )2 


i=-N 


and  minimize  it  with  respect  to  a,  b,  and  c,  where  at  to=0,  a,  b, 
and  c represent  the  estimate  of  position,  velocity,  and  accelera- 
tion respectively.  It  is  then  possible  to  solve  for  the  quanti- 
ties a,  b,  c using  simultaneous  equations. 

There  are  two  advantages  in  using  this  procedure.  First,  the  method 
does  not  require  equally  spaced  data  (frequently,  one  or  more  posi- 
tions in  the  series  is  missing),  although  a test  is  made  for  adequate 
spread  of  the  data.  Second,  by  simply  changing  the  length  of  the 
fit  (2N+1),  we  can  change  the  bandwidth  of  the  low  pass  filter. 

This  is  useful  when  floats  undergo  rapid  velocity  changes  (see 
Study  46). 

The  frequency  response  of  this  operation  is  easily  shown.  Assuming 
no  missing  data,  the  least  squares  fit  gives  the  velocity  estimate 

b=Etlxi/Eti,  where  the  sum  is  from  -N  to  +N.  Since  this  is  an 
antisymmetric  function,  let  x-j=A  sin  toti  where  w is  the  frequency. 

By  writing  u>  = w'  where  2n/2M  is  the  Nyquist  frequency,  and 
ti  = iAt,  we  obtain  b = AE  i s i n ( co ' it i )/AtEi 2 . Thus, 


b*  = £ i sin (oj'-rri  )/Ei 2 becomes  the  attenuation  coefficient  for  the 
velocity  in  dimensional  units  given  by  A and  At  as  a function  of 
frequency. 

The  non-dimensional  frequency  response  of  this  differentiating 
low  pass  filter,  (b*)2,is  shown  in  Figure  6 for  the  cases,  N=6 
and  N=12.  In  computing  the  velocity  stick  plots  in  this  report  we 
used  N=ll,  which  gives  a 1/2  power  cutoff  at  a period  correspond- 
ing to  % 5 1/2  days. 
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TRAJECTORIES  AND  STICK  PLOTS 

Trajectories  and  stick  plots  for  each  float  are  identified 
by  study  number  since  many  of  the  floats  were  launched  more 
than  once.  The  scale  for  the  trajectories  is  two  inches  per 
degree  of  longitude  for  most  studies,  but  studies  37,  39,  41, 
46,  48,  and  49  are  plotted  one  inch  per  degree  of  longitude. 
No  smoothing  is  done  on  the  data,  and  three  positions  per  day 
are  plotted  on  a Mercator  projection.  A symbol  is  plotted 
for  the  first  fix  of  the  day  and  dots  for  the  third  and  fifth 
fixes.  Float  fixes  taken  from  the  ship's  log  at  launch, 
recovery,  and  subsequent  fixes  made  by  passing  over  the  float 
are  labeled  L,  R,  and  S respectively.  The  months  are  labeled 
with  Roman  numerals  on  the  first  occurrence  of  the  month,  and 
days  are  labeled  with  Arabic  numerals  on  the  7th,  14th,  21st, 
and  28th  of  the  month. 

Velocities  for  the  stick  plots  were  computed  using  all  six 
fixes  per  day.  The  data  was  smoothed  using  a filter  of  23 
points,  so  any  studies  with  less  than  four  days  of  data  are 
omitted.  A time  scale  with  month  and  day  labels  is  given 
below  the  plot  and  a speed  scale  of  centimeters/sec.  is  on 
the  side. 
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TEN  DAY  TRAJECTORIES 

These  charts  cover  the  area  27°N-30°N  and  73°W-68°W  for  the 
period  April  1973  through  January  1974.  This  area  and  time 
period  were  chosen  to  show  the  greatest  number  of  floats 
in  a small  area  so  the  relative  motion  of  floats  over  small 
time  periods  could  easily  be  seen.  The  study  number  shows 
the  position  of  the  float  on  the  first  day  of  the  ten  day 
period.  If  the  float  was  launched  during  the  ten  day 
period,  the  study  number  shows  the  launch  position.  The 
position  of  the  float  at  the  end  of  the  ten  day  period, 
or  at  recovery  or  failure,  is  shown  by  an  arrowhead 
indicating  the  general  direction  of  drift.  Table  1 gives 
launch  and  recovery  information  and  is  intended  for  use  in 
conjunction  with  these  charts. 

Positions  are  plotted  six  times  a day.  No  smoothing  or 
filtering  has  been  done  on  these  positions.  Wiggles  on 
the  trajectories  which  occur  at  regular  intervals  are 
due  to  inertial  oscillations,  and  those  which  are  ir- 
regular are  due  to  errors  in  TOA  calculation,  which, 
in  this  area,  are  usually  caused  by  a weak  float  signal. 
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"SPAGHETTI"  PLOT  WITH  BOTTOM  TOPOGRAPHY 

The  "Spaghetti"  plot  is  a composite  of  all  studies  for  the  period 
September  1972  through  June  1976.  Trajectories  are  unsmoothed 
and  one  position  per  day  is  plotted  on  a Mercator  projection. 
Since  there  is  only  one  position  plotted  per  day,  trajectories 
in  areas  such  as  30°N,  70°W  appear  rough  because  of  high  veloci- 
ties and  rapid  changes  in  velocity. 

The  bottom  topography,  superimposed  on  the  spaghetti  plot  is 
reprinted  from  Pratt  (1968).  Contours  are  drawn  at  100  meter 
intervals  with  a few  additional  contours  at  50  meters.  Sea 
mounts  are  represented  with  a Q , areas  less  than  100  meters 
deep  are  shaded  with  a diagonal  line,  and  land  masses  are 
shaded  with  a speckle  effect. 
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